Capillary wrinkling scaling laws of floating elastic thin films under a drop 
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This paper uses dimensional analysis to define the general expression of the pair (V,£@), and 
identifies the dominant combination-parameters of the capillary wrinkling problem, while it also de- 
termines the dominant parameters of different problems relating to its use. The dimensional analysis 
results reveal that, in general, there are no universal scaling laws for capillary wrinkling. Only for a 
small/moderate deformation, it was found that the wrinkling number N is mainly controlled by the 
ratio of bending stiffness and surface tension, while the wrinkling length £ is controlled by the ratio 
of in-plane stiffness and surface tension. Having linear physical relationship in the case of the small 
deformation, simpler scaling laws are proposed for the pair (NV, 2). The universality of the scaling 
laws, which are verified by the dimensional analysis, will give us more confidence. As a natural 
extension, we gave the pair (N, £) a thin film case made of axisymmetric anisotropic materials. By 
using Tanner’s scaling laws, we obtained dynamical scaling laws for a drop radius and the pair 
(N,£), which shows that the pair (N, @) will fade away with time. Finally, we obtained the pair 


(N, £) within the gravity regime. 
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INTRODUCTION 


The deformation patterns of elastic membranes under 
tension is called wrinkling. Wrinkling, which is caused 
by capillary surface tension, is called capillary wrinkling 
(Figure 1). 
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FIG. 1: a. Typical wrinkling of elastic film; b. Typical capil- 
lary wrinkling of elastic film. The problem is to find the de- 
formation wrinkling pattern pair (N, £), where the wrinkling 
number is N and the wrinkling length is £ 


In recent years wrinkling patterns have drawn partic- 
ular attention [1-11], since it can be a useful tool to infer 
material parameters that might otherwise be inaccessible. 
For example, the commonly observed tearing instability 
of an elastic sheet, adhered to a rigid substrate, can also 
be used to characterize the adhesion energy, which is the 
traction forces that are exerted by fibroblasts during cell 
division to the mechanical properties of the membranes 
themselves (such as Young’s modulus and thickness) [5]. 

There has also been considerable interest in under- 
standing fundamental aspects of wrinkled membranes in- 
cluding wrinkling of a compressed elastic film on a viscous 
layer [1], wrinkling mechanics and the geometry of an e- 
lastic sheet under tension [2, 3], the size and number of 
wrinkles [5], the transition from wrinkling to folding [6], 


the analytic analysis of capillary wrinkling of the circu- 
lar elastic membranes [7], wrinkling of pressurized elastic 
shells [8], wrinkling of a charged elastic film on a viscous 
layer [9] and capillary deformations of bendable films [10]. 


Regarding the capillary wrinkling of thin film, a mile- 
stone has been laid by Huang et al. [5] and Vella et 
al. [7]. Huang et al. [5] experimentally derived at the 
length and numbers of capillary wrinkling, while Vella et 
al. [7] theoretically proved the capillary wrinkling length 
expression obtained by Huang et al. [5] for a small defor- 
mation of a circular film. A novel experiment, combining 
both fundamental and applied aspects of the interaction 
between surface tension and elasticity, was presented by 
Huang et al. [5]. In this experiment, a small liquid drop 
was placed onto an elastic membrane that floated in a 
bath of the same liquid. Before adding the drop, the 
membrane was stretched by the surface tension of the 
liquid bath. Once the drop was added, owing to the op- 
posing tension, the contact line of the drop caused radial 
wrinkles. Experimentally, it was found that the wrin- 
kling length £ = 0.031R(42)?, where r is the radius of 
the drop and y is the surface tension coefficient of the 
liquid-gas interface. This pure empirical relationship is 
confirmed by a theoretical justification of Vella et al. [7], 
who used the Foppl-von Karman thin plate theory. 


Although Huang et al. [5] obtained the formula for the 
pair (N, £), where N was wrinkling number, the ques- 
tion remains are these relationships universally applica- 
ble? If the general situation is not, then under what 
circumstances can it be approximated to the universal 
scale law? How can one promote other materials such as 
axisymmetric anisotropic materials? The obtained pair 
(N, £) was the static result, so how is the scale law of the 


dynamics problem established? If the characteristic scale 
of the drop is greater than the capillary scale x71, grav- 
ity begins to dominate the wrinkling process, so what is 
the scale laws beyond the capillary scale «~!? 

The paper first introduce the research topic before pre- 
senting the general expression of the pair (N, £) by using 
dimensional analysis and by determining the controlling 
parameters of the problem,while proposing linear approx- 
imated scaling laws for the pair based on the understand- 
ing of the test’s data, and then extending to the film made 
of composite and laminated materials. For the spread- 
ing of a drop, dynamics of spreading scaling laws are 
proposed based on Tanner’s law. Finally, the paper con- 
cludes with a formulation of wrinkling scaling laws in the 
gravity regime. 


CAPILLARY WRINKLING SCALING LAWS FOR 
SMALL/MODERATE DEFORMATION 


Generally speaking, both the wrinkling number and 
length are functions of the Young modulus EF, film thick- 
ness h, radius of drop r and surface tension y, namely, 
N = F(E,h,v,7,1r) and L = G(E,h,v,7,1r), where the F 
and G are unknown functions. The dimensions of those 
parameters are listed in Table I below. 


TABLE I: Parameters and Dimensions 
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Since the film will be in a bending and stretching s- 
tate by dimensional analysis [12-14], the above formu- 
la can be further expressed into the following forms: 
N = F(D,K,y,r) and £ = G(D,K,y,r), where the 
bending stiffness D = Eh? /[12(1 — v?)] and the in-plane 
stiffness K = Eh/(1 — v2). They can also be expressed 
in dimensionless format, as follows: 


N = Foy/3, | (1) 
‘= afi, J5 (2) 


These are the general scaling relations for capillary wrin- 
kling of the floating thin film. However, both functions F 
and G were not able to be confirmed by the dimensional 
analysis, before other methods must be used. These rela- 
tions reveal that the pair (N, £) will be controlled by two 
dimensionless parameters, namely r,/7/D and y K/y. 
Regarding the capillary wrinkling of elastic film or 
membranes, Huang et al. [5] conducted a good test. 
From the test data we found that the wrinkling number 
N mainly depends on the bending stiffness D; however, 


the wrinkling length £ mainly depends on the in-plane s- 
tiffness K. Therefore, Eqs. (1) and (2) can be simplified 
as follows: 


N = F(r Dp” (3) 
ere (4) 
r y% 


Theoretically, these are the relations that we obtain from 
dimensional analysis. General speaking, the functions F 
and G are not in power forms, which means that capillary 
wrinkling phenomena has no power laws and/or scaling 
laws. Nevertheless, these relations can still provide some 
useful information such as the fact that dimensionless 
parameter r4/Ņy/D is a control variable for the wrinkling 
number N; while the parameter K/y will be a control 
variable for the wrinkling length Z. 

For a small deformation physical laws must be linear; 
therefore, Eqs. (3) and (4) should be a linear function of 
the controlling parameters as shown below 
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where [Y] represents the integer part of real number Y, 
while the constants Cy, Ce and exponents a, 3 can be de- 
termined by either numerical simulations or experiments. 

For a small deformation from the data fitting, Huang 
et al. [5] found C, = 0.033, Ce = 3.62 and a = 1/2, 8 = 
1/2. Hence, the scaling laws are as follows 
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and the number of wrinkling N 


N= [e CAN . (8) 


It is worth pointing out that Eqs. (7) and (8) are univer- 
sal and are valid for all thin flat film that have a small 
deformation. The scaling laws in Eqs. (7) and (8) can be 
used to facilitate measurements of the bending stiffness 
and in-plane stiffness K. If you can firstly get the N, £ 
and radius r, hence, we can calibrate the D and K as 
follows 


3.62 


D= ay VE (9) 
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K = 918.277(%}. (10) 
r 
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These scaling laws have been theoretically confirmed by 
Vella et al. [7], who used von Karman’s circular plate 
theory and hence derived at the similar scaling laws of 
capillary wrinkling. 

With the help of Eqs. (3) and (4), the above scal- 
ing laws of N and £ can be extended to the film, which 
comprises composites materials and laminated materials 
by simply replacing D and K with a corresponding e- 
quivalent or effective bending stiffness Desf and in-plane 
stiffness Keşf; hence, 


y 1/2 
N = |3.62 (« ) ; (12) 
£ we Kept 1/2 
= 0.033( 3 yee (13) 


These expressions provide a pretty accurate estimation 
for those materials without further investigation. This 
is one of advantages of having universal scaling laws 
that are verified by dimensional analysis, which definitely 
gives us ultimate confidence in such an extension. 

For example, if the thin film is made of axial sym- 
metric orthogonal materials, Derg = Dr and Keff = 


1/2 
K,, we have the pair N = [3.02 (raD) | and 
£ = 0.033(K,„/y)!/2, with the equivalent /effective bend- 


ing stiffness D, = E,„h3/[12(1 — v„vo)] and the equiv- 
alent/effective in-plane stiffness K, = E,h/(1 — vroe), 
radius Young modulus EF, and the Poisson ratio vp, and 


the circle Poisson ratio vp. 


CAPILLARY WRINKLING DYNAMICS OF A 
TOTALLY WETTABLE FILM 


Static wrinkling problem has been investigated, how- 
ever, no dynamical wrinkling has been studied yet. A 
drop on a film surface, within a complete wetting regime 
will slowly spread. Typically, the spreading lasts from a 
few hours for ordinary liquids, to several weeks for highly 
viscous fluids such as heavy silicone oils (Figure 2). 
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FIG. 2: Spreading of a drop on a film surface in a total wet- 
ting regime. Tanner’s law: 0p ~ t-3/19_ In this paper, we 
obtained r ~ ¢73/5 


This dynamic process can be expressed in terms of a 
contact angle 0p, which depends on the spreading time t. 


When surfaces are smooth and clean, and for non-volatile 
liquids, Tanner [16] obtained a remarkable universal law: 
Op ~ £73/10, 

The measurements reveal a highly surprising fact, 
namely that the angle @p is completely independent of 
the spreading parameter S = yso — ys — y as long as 
S is positive, that is to say, as long as we are in a total 
wetting regime, where the three coefficients y are the sur- 
face tensions at the solid/air, solid/liquid, and liquid/air 
interfacer, respectively. This is surprising because the 
force F that acts on the system of interest is essentially 
equal to spreading parameter F = yg9—ys_,—y cos Op & 
yso — ysL —y = S. There are two wetting regimes for 
sessile drops. Partial wetting (S < 0): The drop does 
not spread, but instead forms a spherical cap at equilib- 
rium, resting on the substrate with a contact angle 0p. 
A liquid is said to be ”mostly wetting” when 0p < 7/2, 
and ”mostly non-wetting” when 0p > 1/2. Total wet- 
ting (S > 0) : if the parameter S is positive, the liquid 
spreads completely in order to lower its surface (0p = 0). 
The Law of Young-Dupr’e: ycos6p = Yso — Ys, and 
S = (cos @p — 1). 

The precursor film is evidence of the great force F that 
acts on its boundary. The liquid is rapidly drawn towards 
the periphery in the form of a film whose thickness is 
roughly a pancake’s thickness and, which is defined as 
e = 2k! sin(@p/2) as shown in Figure 3 below. 


(b) 


FIG. 3: (a) Liquid drops of increasing size on a sheet of film. 
Gravity causes the largest drops to flatten. (b) Equilibrium 
of the forces (per unit length of the line of contact) act on 
the edge of a puddle. P = (1/2)pge? = —S is the hydrostatic 
pressure. The equilibrium of forces that act on the line of 
contact, y(1 — cosp) = (1/2)pge”, gives the thickness e = 


2x | sin(0p/2), where the capillary length, .~' = \/7/(pg) 


But, behind the film the forces that are involved are 
quite different. Within the drop are forces of traction 
—ysL — ycos@p, whereas within the film (character- 
ized by a zero angle) there are proper forces ysr + y. 
The net force that acts on the drop is then only F = 
y(1 — cos@p) ~ 470%, velocity V = (V*/61)63,, where 
the dimensionless coefficient 15 < l < 20, V* = y/n, 
and viscosity 7. From conservation of the volume Q = 
(7/4)R?0p of the drop, it is easy to obtain the an- 
gle 0p œ~ (Q1/3/V*)3/10¢-3/10 and the radius R(t) ~ 
Q13 (V* /Q1/3) 1/10 41/10, 


Therefore, the radius of the drop is given by 


OM 3754-815. 
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ra etanĝp xk 0h =K (14) 
This scaling laws of the radius of a drop has not been 
recorded in literature. When one substitutes this into 
Eqs(5,6), one has the wrinkling length 


£ = 0.033(5)!/2g7! (%2 )3/54-3/5 
- i a ho —3/5 (15) 
zZ 0.033(=) K (T) t , 
and the wrinkling number 
N= BOGMIVEIRE POA] 
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The dynamics of both length and number are illustrated 
in Figure 4. 


CAPILLARY WRINKLING WITHIN THE 
GRAVITY REGIME 


From Eqs.(14,15), it is interesting to note that both 
£ and N fade away with spreading time, and will stop 
at critical time te = Q1/3/V* = nQ'/3/+. Beyond the 
critical time, the spreading enters the gravity regime. It 
is important to bear in mind that this equation applies 
only when r is less than the capillary length K7! . When 
r >‘, gravity must be taken into account [15]. 

There exists a particular length, denoted «7! = 
\/7¥/(pg), beyond, which gravity becomes important, and 
is referred to as the capillary length. The length «7! is 
generally of the order of a few mm. If one wants to in- 
crease the length «~!, it is necessary to work in a micro- 
gravity environment or, more simply, to replace air with 
a non-miscible liquid whose density is similar to that of 
the original liquid [15]. 

Gravity is negligible for sizes r < «~+. When this con- 
dition is met, it is as though the liquid is in a zero-gravity 
environment and capillary effects dominate. The oppo- 
site case, when r > Kt, is referred to as the ” gravity” 
regime. 

In the critical situation r = «~! and beyond, the scal- 
ing laws Eq.(12) of the wrinkling length £ will be inde- 
pendent of the surface tension y as follows 


K 


l = 0.033( =)", 17 
( F (17) 
and the wrinkling number N, given by 
3.627 | 
= : (18) 
a 


It is clear that both N and £ are constant within the 
gravity regime. 
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FIG. 4: (a) Capillary wrinkling length dynamics & = 
t/ (ae oka (b) Capillary wrinkling number dy- 
namics N* = N/ [pve] 


CONCLUSION 


This paper has attempted to answer all the previous 
questions. In general, there are no universal scaling laws 
for capillary wrinkling. Only in the case of small and 
moderate deformation, can special universal scaling laws 
be formulated. Regarding the bending and in-plane stiff- 
ness, it was found that the wrinkling number N is mainly 
controlled by the ratio of bending stiffness and surface 
tension, and the wrinkling length £ is controlled by the 
ratio of in-plane stiffness and surface tension. By using 
Tanner’s scaling laws, we obtained dynamic scaling laws 
for a drop radius and the pair (N, £), which shows that 
the pair (N, £) will fade away with the time. Finally, the 
pair (N, @) within the gravity regime was also revealed. 

In summary, the highlights of this paper are listed in 
the Table II and III: 


TABLE II: Static capillary wrinkling pair (N, £) 


Dry surface Wet surface 
(Surface tension dominate) | (Gravity dominate) 
N = (3.62 (r\/F) N = [3827] 
` D vpgD 
£ = 0.033r(#)1/2 £ = 0.033(4£)1/2 
7 pg 


TABLE III: Dynamic capillary wrinkling pair (N, £) 


/3\: _ 
N = [3.62( 5)" Va B( 20 3/104 9/10) 


L = 0.033( 4) K! (225 )3/54-9/5 
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